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The evolution of minor planets in the outer part of the asteroid belt is considered. In the frame- 
work of the semi-averaged elliptic restricted three-dimensional three-body model, the boundary of 
regions of the Hill’s stability is found. As was shown in our work (Gerasimov and Solovaya, 1989) 
the Jacobian integral exists. 

The equations of the motion in the rotating coordinate system have the form: 
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The Jacobian integral is 


\{e+v 2 +?)=[v)-c. 


After averaging, the Jacobian integral in the Sun-centered siderical coordinates will have the fol- 
lowing form: 

+ + (3) 
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cos 0^ = cos u> cos lj n + sin u sin lj * cos i , 
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We use a notation similar to the paper Gerasimov and Solovaya (1989), where R is the Sun- 
asteroid distance, p is the Jupiter-asteroid distance, 0 is the angle Jupiter-Sun-asteroid, p is the 
ratio of Jupiter’s mass to the sum of the masses of Sun and Jupiter, m\ and m 2 are stationary 
masses on the rotating axes, ej is the eccentricity of Jupiter’s orbit, l and b are the longitude 
and latitude of an asteroid, lj is the longitude of Jupiter, a,e,i and u> are the osculating elements 
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Figure 2. The regions of stability and instability of the motion projected to the (a,e) plane. 

of an asteroid, and K(k 2 ) is the elliptic i ntegral . When tj = 0 then K(k 2 ) = tt/ 2. In our 
computation \i - 0.0009538. We obtain the boundary of regions of possible stable motion from 
the equation (3) with a fixed value of C and changing elements of fictitious asteroids in the Jimits 
0 < e < 0.4 and 0° < * < 40°. ,. 

For €j = 0 , the value o f C is equal 1.51 938 for the inte rior po int L , for e, = 0.06 2, C — 1.538, 
respectively^ We decided to use the equation (3) for a model, in whic h the heliocentric elongation 
of an asteroid from Jupiter is equal 180°. The regions of the Hill’s stability and instability are 
shown for the three-dimensional case in Figure 1. On the (a,e) plane these regions are plotted in 
Figure 2. 







Asteroids, Comets, Meteors 1991 


567 



Figure 3. The evolution of the perihelion distance q [AU] and the aphelion distance Q [AU] of 
model asteroids within the period of 10 000 years for the first variant of ej and the commencing 

asteroid inclination i — 0°. Asteroids: 0 — 0.00 (the eccentricity), A — 0.04, H 0.08, X — 0.12, 

O — 0.16, + — 0.20, X — 0.24, Z— 0.28, Y — 0.32, X — 0.36, * — 0.40. 


In the outer part of the asteroid belt Hill’s stable asteroid orbits cannot exist with high in- 
clinations. Among real asteroids of the outer asteroidal belt we found only a few asteroids from 
Ephemerides of Minor Planets for 1991 (Batrakov, 1990), which were above the surface represented 
by Figure 1. They are asteroid 153 Hilda (a = 0.766, e = 0.145, i = 7.85°), asteroid 361 Bononia 
(a = 0.757, e = 0.125, i = 12°. 67), asteroid 1212 Francette ( a = 0.760, e = 0.110, i = 7°. 58), and 
asteroid 1529 Oterma (a = 0.769, e = 0.190, i = 9°.03). According to Milani and Nobili (1985) 
and Schubart (1991) asteroids may exist within the outer asteroid belt if there are librations of 
the argument of perihelion. Then their aphelions will always remain far from Jupiter. Numerical 
investigations in the framework of the restricted three-body problem have been developed by many 
authors. They show that asteroids will be ejected by Jupiter in only a few thousand years, with 
the exception of the stable librators. 

We took a dynamic model of the solar system, consisting of six major planets and eleven 
massless fictitious asteroids located at the similar to the four asteroids sis mentioned above. We 
traced the orbital evolution of these asteroids using the numerical integration program with the 
integrator RA15 (Everhart, 1985). The input data a , e and i were taken from the equation (3) 
for the eccentricity of Jupiter ej = 0 — the first variant, and ej = 0.062 — the second variant. 
The starting value of the argument of perihelion was selected from the equation (4) at 0 — 180°. 
The starting epoch of the numerical integration is March 25, 1991 UT. The results obtained by the 
numerical integration for the groups of asteroids with started values of i equal 0° and 10° and for 
the same starting values of a and e are plotted in Figures 3-4. 
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Figure 4. The evolution of the perihelion distance q [AU] and the aphelion distance Q [AU] of 
model asteroids within the period of 10 000 years for the first variant of e 3 and the commencing 

asteroid inclination i = 10°. Asteroids: © — 0.00 (the eccentricity), A — 0.04, H 0.08, X — 

0.12, O — 0.16, + — 0.20, X— 0.24, Z — 0.28, Y — 0.32, W — 0.36, K — 0.40. 


Asteroids with i = 0° are ejected to Saturn after close encounters with Jupiter. Asteroids 
having the inclination ~ 10-20° and the same commencing values of a and e as of the previous 
case are practically all stable. We see, that the inclination of orbits of asteroids from the outer belt 
has the same role in the three-dimensional model as the libration of the argument of perihelion in 
the tw o-dimensional model. When th e inclina tions of asteroids are ~ 10-20°, aphelions of their 
orbits will always remain far from Jupiter. 
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